Origin of the quasi-universality of the graphene minimal conductivity 
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It is a fact that the minimal conductivity ao of most graphene samples is larger than the well- 
established universal value for ideal graphene 4e 2 /irh; in particular, larger by a factor > tt. Despite 
intense theoretical activity, this fundamental issue has eluded an explanation so far. Here we present 
fully atomistic quantum mechanical estimates of the graphene minimal conductivity where electron- 
electron interactions are considered in the framework of density functional theory. We show the 
first conclusive evidence of the dominant role on the minimal conductivity of charged impurities 
over ripples, which have no visible effect. Furthermore, in combination with the logarithmic scaling 
law for diffusive metallic graphene, we ellucidate the origin of the ubiquitously observed minimal 
conductivity in the range 8e 2 /h > a > 4e 2 /h. 
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I. INTRODUCTION 

The ensuing years after the first reported transport 
measurements of isolated graphene flakes have witnessed 
an intense debate on the origin of the quasi- universal 
value of the minimal conductivity ctq. The fact that ex- 
perimental confirmations of the celebrated universal min- 
imal conductivity 4e 2 /7r/i, expected for ideal graphene, 
are the exception^* 2 - rather than the rule and that most 
graphene samples present larger values, typically > 
4e 2 /h, regardless of varying experimental conditions 3 - - — , 
remains one the major fundamental unresolved questions 
in graphene physics. Although this issue has spurred a 
vast amount of theoretical work 9 - " 19 ! 21 " — , to date, the 
answer remains elusive and no theory has been able to 
render a full picture of this jigsaw puzzle. 

Numerical 9,10 as well as analytical work 11 for simple 
models of disorder and no n- interacting electrons indicate 
that, as long as inter valley scattering is avoided, single- 
parameter scaling applies. This means that a beta func- 
tion f3(ao) = dLn(ao) / dLn(L) exists, where L is the size 
of the sample. In particular, /3(ao) > 0, behaving as 
1/(7t<jo) for ao — >• oo. In other words, neutral graphene 
would not present a metal-insulator transition, being its 
conductivity bounded from below by the universal value 
<Jo = 4e 2 /7r/i for ideal graphene and unbounded from 
above (ao — >• oo) as the disorder is increased. This pre- 
diction ultimately relies on the expectation that inter- 
valley scattering is not activated for most common types 
of (long-ranged) disorder and relevant length scales, a 
fact confirmed by some atomistic simulations 12 , but ques- 
tioned by others 13 . 

Alternatively, there have been efforts to consider 
likely sources of disorder such as nearby charged 
impurities^"—, ripples 2 ^, and resonant scatterers^ 2 - in 
the most possible realistic manner. This usually 
comes at the expense of a rigourous quantum me- 
chanical treatment of transport which renders ao size 
independen t 14 i 15 i 21 ~ —. A number of predictions can be 
made in these semiclassical aproaches, having in com- 



mon, e.g., the fact that ao always decreases with increas- 
ing impurity concentration ni mp in addition to not be- 
ing bounded from above as ni mp — ^ 0. This is at odds 
with the theoretical predictions mentioned in the previ- 
ous paragraph and hardly agrees with available experi- 
mental evidence 3 - - — . The inherent limitations of semiclas- 
sical approaches have been recently appreciate d 18 ^ 9 and 
are at the heart of the discrepancy. One cannot deny, 
however, the insight gained on the actual microscopic 
origin of the minimal conductivity from a realistic treat- 
ment of disorder. Combining a full quantum mechanical 
approach to transport with a realistic description of dis- 
order and screening seems to be the only way to resolve 
the controversy and this is our major contribution in this 
work. 

Our findings can be summarized as follows, (i) First, 
our calculations corroborate previous theoretical work 
and agree with experimental findings such as the ob- 
served linear behavior of the conductivity with electron 
density as well as the observed different mobilities for 
electrons and holes, (ii) Second, charged impurities pre- 
vail over ripples in the experimentally relevant range of 
impurity concentrations and distances to the graphene 
flake, increasing the clean-limit conductivity, (v) Finally, 
this increase is limited to a small percentage of the clean- 
limit value for small systems, which, in combination with 
the logarithmic scaling law predicted in the absence of in- 
ter valley scattering, allows us to explain the narrow range 
of values 8e 2 /h > ao > 4e 2 jh of the observed graphene 
minimal conductivity. 

II. HYDROPHENE: A MINIMAL MODEL 

We consider a model for graphene where the C atoms 
are replaced by hydrogenic atoms. The 7r orbit als of the 
C atoms, responsible for the low energy physics, are rep- 
resented by the s orbitals of the hydrogenic atoms, which, 
in turn, are approximated by a gaussian function. The 
gaussian exponent is optimized to reproduce the low- 
energy band structure of real graphene when both are 
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FIG. 1: Color online, (a) A hydrophene ribbon contacted to 
metallic electrodes as to drive the current in the armchair di- 
rection (only the first layer of the electrodes is shown). The 
electrodes are modeled here by a bidimensional three-fold co- 
ordinated tight-binding Bet he lattice, schematically shown in 
the inset, (b) Band structure (in atomic units) in the lo- 
cal density approximation near one of the Dirac points for 
graphene and hydrophene where a minimal basis set has been 
used for graphene. (c) Electrodes density of states per atom 
for a nearest-neighbors hopping amplitude t = 2.7 eV. 



computed in the local density approximation (LDA) as 
implemented in GAUSSIAN03^ [see Fig. [T^b)]. We will 
refer to this model as hydrophene. We are interested in 
the conductivity, defined through a = G^f, where G is 
the conductance of finite hydrophene ribbons of width 
W along the zigzag direction (x) and length L along the 
armchair direction (y). The latter coincides here with the 
direction of the current which is driven by metallic elec- 
trodes contacted along the width of the ribbons as shown 
in Fig. []Ja). The metallic electrodes are modeled by a 
bidimensional tight-binding Bethe lattice of coordination 
three and only nearest-neighbors hopping where the six 
crystallographic directions coincide with those of the rib- 
bon (see inset). This electrode model relates closely to 
graphene, but presents a finite density of states at the 
neutrality point E = [see Fig. [U^c)]. As shown in Fig. 
QJa), a branch of the Bethe lattice is connected to each 
undercoordinated atom on the zigzag edge. The minimal 
conductivity (defined in the limit W/L — » oo) depends 
on the particular choice for the electrode model, but 
never exceeds the universal value an = ^—ItLQ^ for ideal 

u 7rh 

graphene. We have chosen here metallic ribbons, i.e., 
those that, when their electronic structure is computed at 
the simplest nearest-neighbors tight-binding model, this 
does not present a gap for L — >• oo. The results do not de- 
pend either on this choice or the direction of the injected 
current. 




FIG. 2: Color online. Upper panel shows an example of the 
induced charge density due to the ripples shown in the lower 
panel. The scale in the upper panel refers to charge per atom 
which has been smoothed over a fine grid. The scale in the 
lower panel is given in A. 



III. POSSIBLE SOURCES OF DISORDER 

We consider here two types of disorder currently ac- 
cepted to possibly influence the mobility and the min- 
imal conductivity of graphene: Charged impurities and 
ripples. 

The model for the ripples consists of a randomly gen- 
erated position-dependent height function given by the 
equation 



Att 



h(x, y) = V -zr sin (k x x + k y y), 

t N k\k\ 



(1) 



where Nk is the number of Fourier components (~ 5) 
and A is the parameter that accounts for the overall de- 
viation from planarity which is chosen as to reproduce 
typical ripple height-to-size ratios in the order of ~ 0.1. 
2n/W < k x <C 2ix I a and 2tt/L < k y <C 2tt / 'a are random 
reciprocal wave vectors where a is the graphene lattice 
constant. Since in the hydrophene model there is only 
one spherical orbital per site, changes in height only af- 
fect the hopping between atoms (and the overlap), creat- 
ing an accompanying pseudo-magnetic field landscape 26 
which does not affect the charge density distribution. To 
account for changes in the on-site energies due to the sp 3 
rehybridization in real graphene, we have considered a 
scalar potential of the type (/>(#, y) = —B (V 2 h(x 1 y)) — , 
where B > is chosen as to reproduce potential fluctu- 
ations in the order of a few tens of meV 27 . An example 
of the corrugation and LDA induced charge is shown in 

Fig.m 

For simplicity in the analysis of the results, we con- 
sider Z = 1 impurities randomly distributed in the same 
plane at a distance D off the hydrophene ribbon. This is 
probably a good assumption since charged impurities are 




FIG. 3: Color online. Screening electron density induced by 
charged impurities (black dots) placed at a distance D = 10 
(a), D = 5 (b), and D = 1 (c) A. Overall charge neutrality 
within the graphene flake has been imposed in all cases. The 
scale refers to charge per atom which has been smoothed over 
a fine grid. 



expected to be located near the surface of the Si02 sub- 
strate or in between the surface and the graphene flake 
or even adsorbed on the graphene flake, depending on 
their specific origin. Figure [3] shows the screening charge 
density for various representative disorder realizations of 
charged impurities at different distances off the graphene 
plane. Starting at D = 10 A, a landscape of electron- 
hole puddles gives way to strongly localized screening 
charge around each impurity as D — >• 0. This screen- 
ing charge approaches the S— function predicted by effec- 
tive models^, although regularized by the finite value of 
D and by the lattice constant of our atomistic model. 
Figure |4] shows the accompanying Kohn-Sham potential. 
The range of the potential induced by a single impurity 
can be as small as £ ~ 0.5 nm and the depth as high as 
2 eV for D = 1 A. which, in principle, can induce strong 
intervalley scattering 13 . Electron- hole puddles are also 
induced by the ripples due to changes in second-neighbor 
hopping, rehybridization, and local changes in exchange 
interactions 32 , but the average induced electronic density 
is still typically one order of magnitude smaller than the 
one induced by impurities at D = 10 A. 



FIG. 4: Color online. Kohn-Sham potential in the graphene 
plane as induced by charged impurities (black dots) placed at 
a distance D = 10 (a), D = 5 (b), and D = 1 (c) A off the 
graphene plane. The scale is in eV. Overall charge neutrality 
has been imposed in all cases within the graphene flake. 



IV. CONDUCTIVITY RESULTS 

For the conductivity calculations we have employed the 
code ANT.G03, which is part of the quantum transport 
computational toolbox ALACANT— . The basics of 
the calculation are as follows: i) The Bet he lattice is in- 
corporated into the partially pre-computed Green's func- 
tion of the isolated hydrophene ribbon, G#, in the form 
of a self-energy E, ii) a density matrix is obtained from 
the new Green's function imposing overall charge neutral- 
ity in the ribbon, iii) a new Green's function is evaluated 
from the previously obtained density matrix, and iv) the 
procedure is repeated until a self-consistent solution is 
reached. For the LDA exchange correlation functional 
we have employed the standard approximation as im- 
plemented in GAUSSIAN03 20 . With the self-consistent 
Green's function and the self-energies of the left(E^) and 
right (E^) electrodes, the conductance can be calculated 
from the Landauer formalism: 



G(E) 



2e 5 



Tr 



G ] h {E)T r (E)Gh(E)T l {E) 



(2) 



where V 



R(L) 



i E 



R(L) 



E f 



The LDA conductivity a as a function of energy is 
presented in Fig. [5] for a ribbon with an aspect ratio 
W/L « 3. The length here is L w 4 nm and will re- 
main the same in all the calculations. The solid line 
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FIG. 5: Color online. Conductivity as a function of energy. 
The result for a clean ribbon with aspect ratio W/L « 3 is 
shown by the black solid line. The results for various realiza- 
tions of randomly distributed impurities at a distance D = 1 
A corresponding to a density of ni mp ~ 2.0 10 13 cm -2 are 
shown by blue dashed lines. Inset: Same as in main plot, but 
after averaging over disorder realizations. The minima have 
been offset to the origin here. 



corresponds to the ideal or "clean" case with neither rip- 
ples nor impurities. The Dirac point E$ appears shifted 
to positive energies (the Fermi energy has been set to 
zero) due to the influence of the metallic electrodes, de- 
spite the fact that overall charge neutrality is imposed 
on the ribbon. We have checked, by shifting the Fermi 
energy upwards to the Dirac point Eq on electron doping, 
that the overall conductivity curve is not appreciably af- 
fected with respect to the undoped case. From now on 
we will take <Jo = ct(Eq) (for Eq typically > 0) as the 
minimal conductivity. This is now shown as a function 
of W/L (black dots) in Fig. [6] The fact that do scales 
with W/L to a value slightly smaller than 4e 2 /irh for 
W/L —> oc 24 can, in principle, be attributed to the cho- 
sen electrode model. One could possibly improve this 
result, i.e., increase the conductivity closer to the univer- 
sal value 4e 2 /7r/z, by tuning the tight-binding parameters 
of the Bethe lattice or by using a different model for the 
electrode o 24 i 25 . Whether or not the LDA minimal con- 
ductivity of clean graphene can reach the universal value 
4e 2 /7rh is, anyhow, not essential in the ensuing discus- 
sion. 

We now turn our discussion to the effect of charged im- 
purities and ripples on the conductivity. Figure [5] shows 
<j(E) for various realizations of randomly distributed im- 
purities at fixed ni mp and D (dashed lines). Defining 
<j*(E) = a(E+Eo) — (Jo and averaging over impurity real- 
izations, it becomes apparent (see inset) that (a*) oc E a 
with a ~ 1 (more visible for electrons) in contrast to the 
clean limit case where a > 1. This result, in addition 
to our numerical evidence that the chemical potential 
depends linearly on the density for high values of ni mp 
(see also Ref. [33), provides further confirmation that 
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FIG. 6: Color online. Conductivity as a function of the as- 
pect ratio, (a) The results for flat and clean graphene are 
represented by black dots and the ones including only ripples 
by crosses. The results for charged impurities placed at a dis- 
tance D = 5 A off the graphene plane are shown by triangles 
and stars for two different concentrations, (b) Same as in (a) 
but for a distance D = 1 A. 



<j(n) oc n, as experimentally observed and previously ex- 
plained in the Boltzmann transport approximation^"—. 
It also becomes apparent in the inset that the mobility 
(/i = a/ne) of electrons decreases as compared to that 
of holes 34 . All these results are nicely compatible with 
previous works and give us confidence on the validity of 
our LDA results for hydrophene. 

On top of the clean graphene minimal conductivity, 
Fig. [6] also shows (cr ) for a large set of impurity and 
corrugated graphene realizations for two values of D and 
two values of ni mp . Each point has been obtained af- 
ter averaging over 15-20 realizations. The results clearly 
indicate that, for the chosen range of D and ni mp , the 
latter typically one order of magnitude larger than the 
value estimated for exfoliated graphene deposited on 
Si02 (~ 10 11 — 10 12 cm 2 ), charged impurities increase 
the minimal conductivity. Ripples, if anything, add some 
dispersion to the clean-limit values but their influence is 
negligible compared to that of charges, at least in the 
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FIG. 7: Color online. Dependence of the minimal conduc- 
tivity on disorder for a representative ribbon with an aspect 
ratio of W/L « 2.6. (a) Results for two different concentra- 
tions of charged impurities as a function of their distance to 
the graphene plane, (b) Results for two different distances as 
a function of the concentration of impurities. 



range of parameters shown. Importantly, the unavoid- 
able and large statistical uncertainty does not withstand 
the fact that (ctq) scales to a constant value with W/L, 
allowing us to define a minimal conductivity for a given 
length L. A downward deviation is, nevertheless, appar- 
ent at large W/L for strong disorder (D = 1 A) which 
signals the activation of localization. For completeness, 
we also present a systematic study of (ctq) for W/L ^2.6. 
(do) increases as D decreases down to the smallest mean- 
ingful value, D = 1 A [see Fig. Efa,)], and increases 
with increasing ni mp up to the highest value considered 
n imp « 2.5 10 13 cm" 2 [see Fig. [7(b)]. 

A scaling analysis with L is beyond present numer- 
ical capabilities for self-consistent calculations. Never- 
theless, an increasing conductivity with increasing dis- 
order strength is compatible with the scaling law for 
diffusive metallic graphene in the absence of intervalley 
scattering ^ 11 ! 35 : 



a /(4e 2 //i) = -Ln(L/Z intr a), 



(3) 



where /intra can be interpreted as the disorder-dependent 
intra-valley mean free path which, in principle, can be 
estimated from our numerics. Notice that regardless of 
whether or not the weak disorder cases (large D and/or 
small nimp) lie outside the realm of the diffusive regime 
for the system sizes considered, Eq. [3] is certainly valid 
for the strong disorder case close to the critical value of 
W/L beyond which inter- valley scattering decreases the 
conductivity [see Fig. [6(b)]. The key observation is now 
that do is only increased by as much as « 10% with re- 
spect to the clean limit value (see Figs. [6] and [7]) before 
localization sets in. This sets a minimum effective value 
for /intra of ~ 1.5nm with a very weak dependence on dis- 
order over approximately a decade of (smaller) impurity 
concentrations and (larger) distances of the impurities 
to the graphene flake 36 . This minimum length can be 
used now to estimate the maximum possible value of the 
minimal conductivity for the experimentally largest sys- 
tem sizes, (Tq 13 *^ ~ 1/im) « 8e 2 /h. (For larger samples 
L may always be effectively limited by the phase coher- 
ence length < 1/im 37 .) In summary, the weak loga- 
rithmic L- dependence of ao in Eq. [3] and the effective 
minimum intra-valley mean free path in the order of « 
1.5 nm combine to approximately cancel the factor tt for 
relevant length scales (0.1fim< L,l$ < lfim) and explain 
why most graphene samples exhibit minimal conductivity 
values in the range Se 2 /h > <Jo > 4e 2 /h. 

A few final remarks are in order, (i) For strong dis- 
order or very large samples with unintentional disorder, 
i.e., for /inter <C L, where /inter is the inter- valley scatter- 
ing length, graphene behaves as an insulator. The former 
condition may be achieved by intentional doping 8 . The 
latter, however, may be prevented by the temperature- 
dependent , which effectively limits L to < ljam^. 
(ii) When k n tra > L, for instance, for very weak dis- 
order or for very short samples with unintentional disor- 
der, the scaling law in Eq. [3] for diffusive systems does 
no longer apply and <Jo approaches the universal value 
Ae 2 /ith as, e.g., reported in Refs. [l||2|. (hi) Screening from 
the substrate or the occupied bands (not included in the 
model) can only increase /intra, strengthening the quasi- 
universal character of cr . (iv) Finally, when zero-energy 
states (or resonant states) appear due to the presence of 
covalently-bonded adsorbates or vacancies, a number of 
predictions that range from the expectation of an insu- 
lating behaviour— at relevant scales to a finite minimal 
conductivity 2 - 2 have been put forward. These, certainly, 
might play a role in the mobility 39 and minimal conduc- 
tivity of graphene, but this study is out the scope of this 
work. 

In summary, we have evaluated the conductivity of 
graphene in the presence of charge impurities and ripples 
including the screening in the local density approxima- 
tion. Impurities are solely responsible for the increase of 
the minimal conductivity with respect to the clean-limit 
universal value. We have quantified this increase and es- 
timated that the minimal conductivity normally lies in 
the range Se 2 /h > ctq > Ae 2 /h as observed in experi- 
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merits. 
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